Quantum estimation via minimum Kullback entropy principle 
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We address quantum estimation in situations where one has at disposal data from the measurement of an 
incomplete set of observables and some a priori information on the state itself. By expressing the a priori 
information in terms of a bias toward a given state, the problem may be faced by minimizing the quantum relative 
entropy (Kullback entropy) with the constraint of reproducing the data. We exploit the resulting minimum 
Kullback entropy principle for the estimation of a quantum state from the measurement of a single observable, 
either from the sole mean value or from the complete probability distribution, and apply it as a tool for the 
estimation of weak Hamiltonian processes. Qubit and harmonic oscillator systems are analyzed in some details. 

PACS numbers: 03.65.Wj, 03.65.Ta, 02.50.Tt 



I. INTRODUCTION 

Quantum estimation of states and operations is a relevant 
topic in the broad field of quantum information science JTj]. 
The subject has a fundamental interest of its own, since it 
concerns the characterization of the the basic objects of the 
quantum description of physical systems. In addition, quan- 
tum estimation techniques have been receiving attention for 
their role in the characterization of gates and registers at the 
quantum level, which itself is a basic tool in the development 
of quantum information technology. 

In order to characterize a quantum system one may measure 
an observable, or a set of observables, on repeated prepara- 
tions of the system. As a matter of fact, the set of observable 
is usually incomplete, i.e., it is not sufficient to give a complete 
quantum information on the system. In other words, it is not 
possible to deterministically reconstruct the full density ma- 
trix of the system from the measured data |2|. In these cases, 
the question is not that of finding the actual state of the system, 
but rather that of estimating the state that best represents the 
knowledge we have acquired about the system from the mea- 
sured data [3]. If we assume not to have any prior informa- 
tion about the system, and quantify this ignorance by entropy, 
than the best estimate may be found by the Jaynes maximum- 
entropy principle (MaxEnt) [4], which include the information 
obtained by measurements while not allowing one to draw any 
conclusions not warranted by the data themselves. 

On the other hand, in most cases one has some a priori 
information on the state of system under investigation. This 
may come, e.g., by some energy constraint or by the consid- 
eration that the system has been weakly perturbed from an 
initial given preparation which is under control of the experi- 
menter. A question naturally arises on how this a priori infor- 
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mation may be incorporated into the estimation procedure and 
whether this can be done together with constraint of reproduc- 
ing the observed data. The quantum mechanical way to incor- 
porate some a priori information is that of bias of the state 
to be estimated towards a given quantum state r, which ex- 
presses the information at hand of the experimenter and which 
may contain one or more free parameters depending on the 
amount of a priori information at disposal. Upon quantifying 
the bias towards r through relative (Kullback) entropy |5l |6[] 
an estimate may be found by minimizing this quantity with the 
constraint of reproducing the data. This minimum Kullback- 
entropy principle (mKE) emerges naturally as a way to es- 
timate the quantum state of a system from incomplete data, 
when some a priori information on the state is available, i.e., 
when a bias towards a known quantum state is present. We 
will exploit this ideas to estimate the quantum state of a sys- 
tem from the measurement of a single observable, either from 
the sole mean value or from the full probability distribution, 
and apply it to the estimation of weak Hamiltonian processes. 
In particular, the case of a weak, but otherwise generic, Hamil- 
tonian is analyzed in some details, with emphasis to qubit and 
harmonic oscillator systems. 



The paper is structured as follows. In Section|II]we review 
classical Kullback-Leibler divergence for probability distribu- 
tion and its quantum counterpart, the quantum relative Kull- 
back entropy. In addition, we state the minimum Kullback 
entropy principle which is used in the subsequent Sections. In 
Section|lII]we exploit mKE for state estimation. In particular, 
we analyze mKE estimation from the measurement of a sin- 
gle observable, either from the sole mean value or from the 
complete probability distribution. Qubit and harmonic oscil- 
lator systems are analyzed in some details. In Section HVl we 
address mKE as a tool in order to estimate a weak Hamilto- 
nian through suitable measurements performed on the evolved 
states. Section [V] closes the paper with some concluding re- 
marks. 
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H. KULLBACK-LEIBLER DIVERGENCE AND QUANTUM 
KULLBACK ENTROPY 



Let us consider a classical system which can be in any of N 
states and let p k be the probability for the k-th state. Suppose 
we want to estimate the probability distribution p = {p k }. If 
all that we know is the number of possible states, there is no 
way to do better than choosing p k = l/N, Vfc. Fortunately, 
in general we have further information about our system, e.g., 
we have at disposal the value of certain functions of the state 
of the system. Thus, we can use these as constraints to es- 
timate p. The problem of finding the most likely distribution 
satisfying a set of given constraints was solved by Jaynes, who 
proposed a general method of inference known as principle of 
maximum entropy (MaxEnt) J3l ■ According to MaxEnt prin- 
ciple the best estimate of p is that maximizing the Shannon 
entropy 



H(p) = -^p fe logp fe 



(1) 



given the constraints. As it was described by Jaynes, such a 
distribution is the one that agrees with what is known, but ex- 
presses maximum uncertainty with respect to all other matters 

a. 

On the other hand, it is not unlikely that the probability 
distribution we want to estimate is biased towards a prior dis- 
tribution q = {q k }- In this case it is useful to introduce the 
Kullback entropy Jill]: 



K(p\q) = ^Pk log{pk/q k ), 



(2) 



that represents the relative entropy of a posterior distribution 
p relative to the prior q. It is worth noting that in the absence 
of any knowledge on q, i.e., qk — l/N (uniform distribution), 
we have K (p\q) = log N — H(p): this is just the opposite of 
the Shannon entropy up to an additive constant. On the con- 
trary, when q is accessible, the Kullback entropy (f2]i leads to a 
generalized method of inference: the most likely distribution 
p with a bias towards q is the one minimizing K(p\q) given 
the constraints. This is the classical principle of minimum 
Kullback entropy (mKE), also referred to as the principle of 
minimum relative entropy ifsL f^L fTotl -, w hich found applications 
in several branches of science Illlll2l[l3ll . Indeed, minimiz- 
ing the relative entropy has all the important attributes of the 
maximum entropy approach with the advantage that prior in- 
formation may be easily included. 

In the following we will deal with the quantum version of 
the mKE principle and its application to quantum estimation 
of states and operations. Before going to the main point let us 
spend few words on how this works in the classical case. Let 
us consider the problem of finding the posterior p given the 
prior q and the constraints J2 k Pk = 1 (normalization) and 
J^kPkAk — (A), i.e., we assume to know the first moment 
of the quantity A. Using the Lagrange multipliers method the 



problem reduces to the minimization of 



F(p,\,r,) =K(p\q 




+ A I ^2,PkAk - (A) 



(3) 



T] and A being Lagrange multipliers. We have: 

= \og( P k/qk) + l + A k \ + r 1 (k = l,...,N) (4) 

1 = 5>fc (5) 



(A) = ^2 Pk A k . 



(6) 



Solving the first equation with respect to p k and using the con- 
dition (0 we obtain: 



Pfc(A) 



<7fce 



-A k x 



-A,\ 



(7) 



where we explicitly wrote the dependence of p k on the La- 
grange multiplier A. Note that for A = one has p k = q k - 
Differentiating Eq. (0 and using the constraint © we obtain 
the differential equation: 



dp k 
dX 



-(A k -(A))p k . 



(8) 



Upon considering the distribution p as a point in the probabil- 
ity distribution simplex, then Eq. ^ can be seen as a trajec- 
tory in such a space. According to the Lagrange multipliers 
method, if we integrate the trajectory (O assuming p(0) = q 
for A = 0, then the minimum of Kullback entropy is achieved 
when the trajectory passes through the surface satisfying the 
constraint EfcPfc^fc = ( A ) & 



A. Minimum Kullback entropy principle 

Let us turn our attention to quantum mechanics. The prob- 
lem is now to find the most appropriate density matrix g de- 
scribing a quantum system that should satisfy some given con- 
straints, which, in turn express the results of an incomplete set 
of measurements. If there is no a priori information, then the 
optimal choice is given by the density matrix g which maxi- 
mizes the von Neumann entropy 



H(g) = -Ti[g\ogg], 



(9) 



and satisfies the constraints, i.e., reproduces the observed data. 
This is the quantum counterpart of the MaxEnt principle intro- 
duced above [3]. On the other hand, when we have some a pri- 
ori information about the system under investigation then the 
state to be estimated has a bias towards a prior one r. This is 
the case, for example, of a quantum system evolved according 
to an unknown but weak Hamiltonian starting from a known 
initial state. In order to build a proper estimation scheme for 
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these situations let us consider the quantum Kullback entropy, 
defined as follows OUll]: 



K{q\t) = Tr[ (loge-logr)] 



(10) 



that is the relative quantum entropy of the density matrix g 
with respect to r. 

As well as in the classical case a probability distribution 
can be seen as a point, in the quantum case the role of "point" 
is played by the density matrix g. We can also associate a 
family of surfaces Tr[Ag] = (A) to each quantum observable 
A. Moreover, it is still possible to define a suitable metric in 
the space of the density matrices (the Hilbert space) lfl5ll : this 
allows us to introduce the infinitesimal increment dg. In this 
way one can obtain the following equation for the quantum 
trajectory Q: 



(ii) 



where {X, Y} = XY + YX and A is again a Lagrange mul- 
tiplier. Eq. (fTTT i is the quantum counterpart of Eq. ([8]) and sat- 
isfies the constraints Tr[g] = 1 and Tr[Ag] = (A), i.e., we are 
assuming that the expectation value (A) of A is known. By 
formal integration of Eq. (TTTT i, it is straightforward to obtain 
the solution 



0(A) = 



,-AX/2 T e -AX/2 
Tr[re~ AA ] : 



(12) 



where we assumed that the formal integration starts from 
g(0) = t, t being the prior density matrix. 

There are cases in which the trajectory ( fTTT ) yields to the 
optimal state with respect to the quantum Kullback entropy, 
i.e., K(q(X)\t) is minimized IU5I1 . In this paper we focus our 
attention precisely on one of these cases, namely, when the 
prior t and the posterior g are close each other according to 
the Fisher information metric. 



B. Remarks 

As it may be expected, when the prior information is very 
weak mKE principle reduces to MaxEnt one. On the other 
hand, some prior knowledge is often present and mKE prin- 
ciple fully exploit the additional information to improve the 
estimation procedure. A question may arise on the choice of 
the measurements: as a general rule those should be tuned in 
order to add information with respect to the prior being other- 
wise useless for estimation purposes. In other words, the in- 
formation coming from the measurements should not be sub- 
sumed by the prior one. 

Kullback-Leibler divergence and quantum Kullback en- 
tropy are involved in different aspects of quantum estimation 
of states and operations, including assessments of priors [ 16] 
in the context of quantum Bayes rule B17I1 . Here we briefly 
mention few relevant applications in order to enlighten analo- 
gies and differences with our approach. 

Let us first consider a situation in which we have no a pri- 
ori information and want to estimate the state of a system from 



the measurement of a set of a set of projector Ak = \<fk) {fk\- 
In this case the maximum-likelihood estimate of the state 111811 
is the density matrix g that minimizes the Kullback-Leibler 
divergence ^ fe pf. log(pfe/ gk) between the observed probabil- 
ity distribution and the quantum mechanical prediction gk = 
(<fk\Q\fk) EH. In other words, maximum-likelihood estima- 
tion leads to the state that fits the given data obtained by the 
given measurement without using prior information about the 
quantum state. 

In Ref. 12011 Kullback entropy is used for quantum esti- 
mation as a loss (cost) function in the search of the optimal 
predictive density matrix in a Bayesian (global) approach. In 
other words, best estimate is found by minimizing the average 
Kullback entropy with respect to the true state. In the same 
perspective, Kullback entropy has been also used as a regu- 
larizing functional in seeking solutions to multivariable and 
multidimensional moment problems l2lll . 

Finally, notice that symmetrized version of the Kullback 
entropy has been also suggested 12211 and may be employed to 
assess the distance between quantum states. 



III. MINIMUM KULLBACK PRINCIPLE FOR STATE 
ESTIMATION 

In this section we exploit mKE principle for the estimation 
of the full density matrix g from an incomplete set of mea- 
surements. We assume that the system under investigation 
has a bias towards the known prior t and first consider to have 
access to the mean value of a single observable. Then the 
analysis is generalized to the case of N observables, with ap- 
plication to the measurement of the complete distribution of a 
single observable. 



A. Measurement of a single observable 

In this Section we assume that only the mean value of the 
observable A can be measured. This is the simplest obser- 
vation level [23] one may devise for a quantum system, and 
is generally considered to provide only a little amount of in- 
formation about the state under investigation. Upon applying 
the quantum mKE principle, the best estimate for the density 
matrix, compatible with the bias is given by Eq. (fl2b . By in- 
troducing the partition function Z = Tr[re~ ] the estimated 
density matrix reads as follows 



z 



(13) 



Moreover, using the spectral decomposition A = 
ak\tpk){<Pk\> {l^fc)} being a complete orthonormal 
system of eigenvectors of the operator A, we can write: 
exp{-iAA} = J2k ex P{ _ 5 a fcA} so that 



(<p m \Q\V n ) = V^ +Q ") A (^ m |T|<p n > 



(14) 



(15) 
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In this way, given the initial density matrix r it is possible to 
estimate the complete state g as follows 



Q : 



<Pm){<Pn 



(16) 



where the value of the Lagrange multiplier A is obtained as a 
(numerical) solution of the equation Tr [qA] = (A), i.e., 



(A) = 



EkvV*l r l*>*) e 



-a k X 



Ok 



T, k ( l Pk\T\<Pk}e- 



a k X 



(17) 



which, of course, is equivalent to the relation —d\ log Z = 
(A). It is worth noting that the estimate (IT6b for the den- 
sity matrix has support in the same subspace of Hilbert space 
where r does. Moreover, if r = |<^n)(<^n|, i.e., if the prior 
density matrix is a projector onto the subspace generated by 
the eigenvector \ipn) of A, then (<p n \T\(p n ) — 5 n n and the 
mKE principle reduces to the MaxEnt one, as in the classical 
case. 



B. Measurement of N observables 



Finally, taking matrix elements of (fT9b and using (l22l and (l23l 
it is possible to reconstruct the posterior state, given the initial 
density matrix r and the measured probabilities p k 



n.m 



(<Pm\T~\<p n ) 



in n 1 1~\ fin ( 



\/PmPn \ip m )(lfn\- (24) 



Notice that in this case we have been able to back-substituting 
the Lagrange multipliers, i.e., Eq. d24b no longer depends on 
the A/j's. 



C. The qubit case 

Here we address the estimation of a qubit state starting from 
the measurement of a single observable [26]. In order to ap- 
ply the mKE principle we assume to have a bias towards the 
state r and choose an observable to measure on the system. 
The measured quantity is the spin along direction ft, which is 
described by the operator: 



A = 



(25) 



Here we assume that the mean values of N different observ- 
ables Ah, k = 1, . . . ,N are experimentally accessible. This 
means that we we have N constraints (A k ), to be considered. 
Upon writing the trajectory (fTTT i for each constraint, we have 
a system of N differential equations which can be written in 
the following compact form: 



N 

£ 

fc=i 



dg 



|{ft£(4k-<4b>)}, (18) 



The solution can be written as 

Q = 



(19) 



where Ai, . . . , Xn are Lagrange multipliers, and we assumed 
g(0) = t. The partition function Z reads: 



Z = Tr 



■T, k A k X k 



(20) 



Again, the values of the multipliers A^ are obtained solving 
the system of equations 



Tr[gA k ] = (A k ) 



1. 



, N 



(21) 



The above analysis allows us to apply mKE principle also 
starting from the measurement of the full distributions of an 
observable. In this case the set of observables to be taken 
into account are the orthogonal (commuting) eigenprojectors 
Ak = \ fk){fk\, (fklfs) — dks of the measured observables. 
The constraints Tr[p^4fe] = p k , correspond to the measured 
distribution. Eq. d20i > and (fJTJ rewrite as 



Z 



{(Pk\r\ipk) 



(22) 
(23) 



where we defined the vector er = (tT%, a%, as), cr k , k = 
1,2,3, being the Pauli matrices. Upon writing the prior state 
in the Pauli basis 



r=i(l + f-«r), |f|<l, 



Eq. ( fT3l reads: 



g= -(1 + v-a), 



(26) 



(27) 



with 



t + 2 sinh (A/2) (n • t) n — sinh A n 

v = r\ — - - ■ u x ' *^ 28 ) 

cosn A — t • n smn A 

where we used Z = cosh X — f-n sinh A. Now, thanks to the 
constraint: 



Ti[gri ■ er] 



n • cr 



(29) 



we can calculate the value of the Lagrange multiplier A, ob- 
taining: 



A = arctanh 



1 — (n ■ er)(r • ft) 



(30) 



In order to have a compact expression of the result let us con- 
sider an operator basis composed by spin operators along three 
orthogonal directions Hi _L n 2 -L n 3 , with fti = ft. In this 
way we can express the components of the vector v as follows: 



v • n\ = \n ■ a) 



v ■ n k = t ■ n k 



(71 ■ <T 



l-(f-n) 2 



(fc = 2,3). 



(31) 



(32) 
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Eqs. OTI ) says that the estimated Bloch component in the di- 
rection of the measured observable is equal to the measured 
mean value, whereas the two other orthogonal components are 
obtained from the prior one by a common shrinking factor. 

As an example, let us assume n — (1,0,0) and f = 
(0, 0, 1): this is the case of the measurement of <j\ (spin 
along x direction) and bias towards +z direction. We have 
x = — (<Ti) and, then, 



which satisfies Eqs. OTI ) and d32t . 



(33) 



D. The harmonic oscillator case 

Now we face the problem of estimating the state of a har- 
monic oscillator with a bias towards a coherent state r = 
|a)(a|, with a £ C. Let us consider a photon number mea- 

7tl = 



surement, i.e., if a denotes the annihilation operator, 



1, the observable expressed by A = era 
{\n)} being the photon number basis. Since 



(n\T\m) 



V n\m). 



T,n n \ n )( n l 



(34) 



we have that the mKE estimated state is still a coherent state 
with amplitude 



-A/2 



Now, using the constraint 



Tr[£»a T a] = N, 



(35) 



(36) 



with N obtained from the experiment, we can evaluate the 
value of the Lagrange multiplier A, namely, 



A = log(H 2 /iV) 



(37) 



also obtaining Z = cxp{N}. Finally, upon substituting in 
([Tol l we arrive at 



q = e 



-N 



, 2 \(n+m)/2 a a" 



|ViVe^)(v / 7Ve^|, 



(38) 



with (f> — arg a. In other words, the best estimate according 
to mKE is a coherent state with average number of photons 
equal to the measured one and phase equal to that of the prior 
coherent state. Notice that the best estimate obtained using 
the MaxEnt principle with the same constraint on the average 
number of photons, but without the bias, would have been a 
thermal state with N thermal photons. 

If by some means the complete photon distribution p n is 
available, the reconstructed state, given by Eq. d24l . reads as 
follows (we assumed that the bias is still towards r = |a)(o;|, 
a e <D): 



Pn 



n)(m\ 



(39) 



Remarkably, Eq. d39l no longer depends on the amplitude of 
the prior state, which enters only trough the phase <j>. This 
makes the above scheme quite promising though measuring 
the photon distribution is, in general, a challenging task. On 
the other hand, in the optical case it is possible to reconstruct 
the p n by means of on/off photodetection and maximum like- 
lihood algorithm l27l l28l]. a method that has been recently 
verified in laboratory 11291 13011 . Multichannel fiber loop detec- 
tor i3~ll l32ll may be also used. 

Finally, we notice that a special case of bias is that toward 
a Gaussian state. In fact the Kullback relative entropy of a 
state q with respect to a Gaussian state r, with the same co- 
variance matrix reduces to the difference of the Von Neumann 
entropies lEBTl 

K(g\r) = H(t) - H{q) 

and thus the mKE principle reduces to the MaxEnt. Notice, 
however, that this is not in contrast with the results above, 
since in that case the results of the measurement do not impose 
the equality of the covariance matrices. 



IV. MKE ESTIMATION OF WEAK HAMILTONIANS 

In the previous section we have shown how to fully recon- 
struct the density matrix of a quantum system from incomplete 
data. Here we will see how it is possible to estimate a weak 
Hamiltonian H by means of the mKE and suitable measure- 
ments onto the evolved states. The idea behind this method 
is that, upon considering weak Hamiltonian processes, the 
evolved state is not too different from the initial one, i.e., it has 
a natural bias towards the unperturbed state. This allows one 
to estimate the parameters (matrix elements) of the Hamilto- 
nian from data obtained by an incomplete set of measurements 
on the evolved state, i.e., to use mKE principle as a effective 
tool for process estimation. 



A. The qubit case 

The initial state r and the evolved state of the qubit systems 
under investigation are connected by the transformation 



et =e- im re im 



(40) 

In. Eq. d40l t is the time evolution. Using the Pauli basis, we 
can express the initial state and the Hamiltonian as 



and 



T=i(l + f-<x), |f|<l, 



ff = ^ri„<7„, |h| < 1, 

u=0 



(41) 



(42) 



respectively, where <x„ are the Pauli matrices with er = 1, 
and er has been defined above. Expanding at the first order in 
h Eq. ( |401 l, we obtain: 



Q t =T + U[T, H]+o{\h\ 2 ), 



(43) 
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with 



[t, H] = i ^ £ksi cru 



(44) 



k,s=\ 



Eksi being the totally antisymmetric tensor, e 12 3 = 1. In this 
way, expansion (l43l can be written as follows: 



g t = - (1 +w- a), 
where the Bloch vector is given by 

w = f + 2h x f . 



(45) 



(46) 



Eq. d46l ) represent a system of equation for the unknowns h. 
The transfer matrix is singular, but the system may be anyway 
inverted using the Moore-Penrose generalized inverse B24I1 . 
thus leading to the expression 



T X W 



(47) 



which provides an estimate for the Hamiltonian h once an es- 
timate for w is given. The latter is obtained upon the mea- 
surement of a spin observable described by the operator n ■ <x 
onto evolved state g t . The best estimate for w according 
to mKE principle is given by Eq. d28T l. By substituting in 
Eq. ( |47l ) and using ( f30b an estimate for h is achieved: 



1 (1 - Vl - k 2 )t-u- 



2f 



?|2 



where 



1 — t ■ n k 



1 — (ft ■ <?)(t • n) 



(48) 



(49) 



As it is apparent from Eq. (08]) the mKE estimate for the 
Hamiltonian Bloch vector is orthogonal to both the prior state 
one and the measurement direction. By a repeated randomized 
choice of the measurement an effective reconstruction may be 
achieved for any, but weak, qubit Hamiltonian. Generaliza- 
tions to higher dimensional systems may be also designed. 



B. The harmonic oscillator case 

Let us now assume that the expansion 

g(t) =T + it[T,H]+o(\H\ 2 ), (50) 

refers to the state of a harmonic oscillator evolving under the 
action of a weak Hamiltonian H that we want to estimate. We 
also assume that the full distribution of a single observable 
can be measured of the evolved state. Therefore, the evolved 
density matrix g t may be reconstructed by mKE starting from 
the observation level Ak — \<fik)(fk\, as described in Section 



IIII Bl and thus obtaining the state (124-b . Using the basis {| 
we can write Eq. ( |50l as follows: 



ms 1 sn 



(51) 



where we used H mn = (ip m \H\ip n ) . Using the mKE estimate 
(l24l l for the evolved density matrix we obtain the following 
hierarchy of equations 



T 7~m.n. 1 



t 



PnPr, 



{Wn\T\Wn){V m\T\ (Pm ) 



(52) 



where p n are the measured probabilities (the constraints used 
for the mKE) and the matrix elements H nm are the unknowns. 

A relevant example, in which the Hamiltonian can be ef- 
fectively estimated using mKE, is that corresponding to H = 
(ga + h.c), a being the annihilation operator starting from 
the sole measurement of the photon distributions. The evolu- 
tion imposed by the Hamiltonian H corresponds to the uni- 
tary displacement operator D(j3) = exp(/3a^ — /3*a), (3 = gt. 
The problem is then to estimate the displacement amplitude 
(3 from the measured photon distribution. We assume that the 
initial state is a coherent state g(a) = \a) (a\. For the sake of 
simplicity we take a and j3 as real. Using the photon number 
basis the evolved state may be written as 

g(P) = e^ 2 J2 ia+ ^ m \n)(m\. (53) 
V n\ ml 

n,m v 

Assuming that the measurement of the photon number is made 
on the evolved state and that mKE principle is used to esti- 
mate the density matrix we equate the above expression to 
that given in Eq. d39l for tfi — (recall that a and (3 are taken 
as real). We thus obtain the following set of equations: 

- (a + (3f + (n + m) ln(a + (3) = In ^Jn\m\p n p mi (54) 

to be solved for (3. It is worth noting that in order to esti- 
mate (3 one can choose to measure a. finite number of i.e., 
k = 0, . . . , N — 1. As a matter of fact, this choice also se- 
lect a subspace of the Hilbert space where the reconstructed 
state is defined. In turn, d54l i corresponds to N 2 determi- 
nations of the same parameter (3. Notice that without using 
the mKE principle the only way to exploit the information 
at disposal, i.e., the elements of the probability distribution 
p n = e -(" +/3)2 (a + (3) 2n /n\, n = 0, . . . , TAT - 1, is to in- 
vert those relations. In order to estimate (3 one should solve 
the set of equations 

- (a + (3) 2 + 2n ln(a + (3) = ln( n\p n ), (55) 

which provide only N determinations of (3. 



V. CONCLUSIONS 

In this paper we have considered quantum estimation of 
states and weak Hamiltonian operations in situations where 
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one has at disposal data from the measurement of an incom- 
plete set of observables and, at the same time, some a priori 
information on the state itself. By expressing the a priori in- 
formation in terms of a bias toward a given state the best es- 
timate is obtained using the principle of minimum Kullback 
entropy, i.e., by taking the state that reproduces the data while 
minimizing relative entropy with respect to the bias. The 
mKE principle has been used to estimate the quantum state 
from the measurement of a single observable, either from the 
sole mean value or from the complete probability distribution. 
In particular, we have analyzed qubit and harmonic systems 
with some details. We have also considered the problem of 
estimating a weak Hamiltonian processes. In this case there 
is natural bias of the evolved state towards the initial state and 
the mKE principle can be used as a tool to estimate the Hamil- 



tonian from an incomplete set of measurements. 

Overall the minimum Kullback entropy principle appears 
to be a convenient approach for quantum estimation in real- 
istic situations and a useful tool for the estimation of weak 
Hamiltonian processes. 
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